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^' Abstract 

A new fully quantum method describing penetration of packet from internal well outside with its tunneling 
through the barrier of arbitrary shape used in problems of quantum cosmology, is presented. The method 

^^v , allows to determine amplitudes of wave function, penetrability Tbar and reflection -Rbar relatively the barrier 

rvj ' (accuracy of the method: |Tbar + ^bar — 1| < 1 ■ 10~^^), coefficient of penetration (i. e. probability of the 

packet to penetrate from the internal well outside with its tunneling), coefficient of oscillations (describing 
oscillating behavior of the packet inside the internal well). Using the method, evolution of universe in the 
closed Friedmann-Robertson- Walker model with quantization in presence of positive cosmological constant, 
radiation and component of generalize Chaplygin gas is studied. It is established (for the first time): (1) 
oscillating dependence of the penetrability on localization of start of the packet; (2) presence of resonant 
values of energy of radiation -Erad , at which the coefficient of penetration increases strongly. From analysis of 
these results it follows: (1) necessity to introduce initial condition into both non-stationary, and stationary 
quantum models; (2) presence of some definite values for the scale factor a, where start of expansion of 
^ ' universe is the most probable; (3) during expansion of universe in the initial stage its radius is changed not 

r^ , continuously, but passes consequently through definite discrete values and tends to continuous spectrum in 

'^ ' latter time. 

■^ 
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1 Introduction 

" T-^ I 

j^ ' If we analyzed existed variety of quantum models which describe formation of the universe and its subsequent 
JH \ evolution in the first stage, we should come to conclusion that the semiclassical approach for description of 
tunneling and determination of wave function is the most prevailing today. This approach forms a basis, 
which props up both models with the Feynman formalism of path integrals in multidimensional space-time, 
developed by the Cambridge group and other researchers, called the "Hartle-Hawking method" (for example, 
see Ref. [T]), and methods based on direct consideration of tunneling in 4-dimensional Euclidian space-time 
called the "Vilenkin method" (for example, see Refs. [21-[TU]). The models in 4-dimensional space-time directed 
on description of inflation, on study of fluctuations of vacuum with inclusion of massive fields (for example, 
see Refs. [llj), multidimensional cosmological models (for example, see Refs. [12l[T3]), variety of string models 
(for example, see Refs. [IHIIS]) have mainly such a semiclassical grounds. To date, this basis is supposed to 
be sufficiently reliable and give interested quantum characteristics of the universe with good accuracy. Such a 
point of view so has been taken root and is prevailing, that, in spite of almost 40 years of researches in quantum 
cosmology (see the first papers |161I17) ). the papers devoted on more tiny and deeper study of possible quantum 
nature of the universe, its formation and evolution on the first stage, can be found enough rarely (for example, 
see [HHin], also [2122,). 

However, one should ask whether the penetrability determined according to the semiclassical theory by a 
shape of the barrier solely between two turning points, gives exhaustive answer and the best estimations of rates 
of evolution of universe. 

(1) If, despite such widespread confidence in the semiclassical approach, we still want to check it, we imme- 
diately will miss some of the parameters. For example, it would seem, one can use a test of T + _R = 1 (where 
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T and R are penetrability and reflection relatively the barrier in the cosmological problem), where usually we 
have no doubt. However, one should recall that in quantum mechanics, the semiclassical approximation neglects 
the reflected waves completely (see [23], eq. (46.10), p. 205, also p. 221-222) and, therefore, to compare the 
calculated penetrability T is just not with anything. 

(2) If we still wanted to determine the reflection coefflcient, then we should increase the order of approxima- 
tion of the semiclassical method (in order to take into account the decreasing component of the wave function 
on the background of increasing one in the region of tunneling), and just stumble on the next problem — pres- 
ence of a non-zero interference between the incident and reflected waves in the radial task. Now the criterion 
T + R — \ for testing is not satisfied and it needs to take into account the third component M of interference in 
addition (see [24] ) . In particular, at unsuccessfully chosen separation of the exactly known full wave function on 
the incident and reflected waves (but the semiclassical approaches have no needed apparatus for such analysis), 
the interference component can increase without limit, and be substantially larger than the penetrability and 
reflection. After the appearance of such an arbitrariness, the penetrability and reflection can freely exceed the 
unit and increase without limit. In what is now the general meaning of the penetrability? 

(3) We shall give only some easy examples from quantum mechanics, (i) If we consider two-dimensional 
penetration of the packet through the simplest rectangle barrier (with finite size), we shall see that the pene- 
trability is directly dependent on direction of tunneling of the packet. So, the penetrability is not a single value 
but the function, (ii) If to consider one-dimensional tunneling of the packet through the simplest rectangular 
barrier, one can obtain "interference picture" of its amplitude in the transmitted region, which is depended on 
time and space coordinates and is an exact analytical solution. Of course, the stationary part of such a result 
coincides exactly with well known stationary solutions |25| . From the arguments above the impression could 
appear that the penetrability defined solely by the shape of the barrier between two turning points is nothing 
more than prevailing simplified understanding, while for more accurate and deep analysis we need in the strong 
basis. 

(4) Advance of the semiclassical approach is in simplicity of formula of the penetrability based on determi- 
nation of the outgoing wave in the asymptotic region. A tunneling boundary condition [7l[8] seems to be natural 
and clear, where the wave function should represent an outgoing wave at large scale factor a. However, whether 
is such a wave free? In contrast to problems of quantum atomic and nuclear physics, in cosmology we deals 
with potentials, which modules only increase with increasing the scale factor a (also their gradients increase, 
which have sense of force acting on the wave) and, therefore, we have nothing mutual with a free propagation 
of the wave in the asymptotic region. Now it is unclear to which combination two Airy functions should be 
combined at turning point, in order to obtain the proper outgoing wave. It turns out that instead of the free 
wave in the asymptotic (missing in problems of quantum cosmology) , we should be able to work with the waves 
propagating inside strong fields (see [24]). 

These problems violate (destroy) the basis of the semiclassical models, and now statements about reliability 
of the semiclassical models are transformed into the question of " faith" in them, though widespread !24; . The 
semiclassical approach could be compared to "black box" in which deeper and more detailed information about 
the dynamics of the universe is hide. More importantly, in such a black box those missing elements are hidden, 
without which it is impossible to combine everything together to obtain self-consistent formalism of quantum 
description of the formation of the universe and its evolution in the first stage. In order to clarify these questions, 
we have developed a new fully quantum method presented in this paper. 

This paper is organized so. In Sec. 2 a new fully quantum method for description of the formation of the 
universe and its evolution in the first stage on the basis of a packet, which penetrates from the internal potential 
well outside by tunneling through the barrier, is presented. A main advance of this method is determination of 
characteristics of the packet with high accuracy (without implication of the semiclassical approximations) . The 
formalism of the method has been developed relatively the barrier of arbitrary shape, that makes the method 
universal. In Sec. 3 the method is applied for solution of the problem of evolution of the packet in FRW-model 
with radiation and Chaplygin gas. The penetrability and refiection relatively the barrier are calculated. We 
propose new characteristics, more adequately determining the probability of formation of the universe and its 
subsequent expansion. An accuracy of the method is demonstrated, achieving to iTbar + -Rbar — 1| < 1 • 10~^^ 
inside whole under-barrier range of the energy of radiation (M = 0, author has not yet found competitive 
approaches, achieving such a precision), stability in calculations for the obtained results is shown. A special 
attention in analysis is devoted to study of the initial conditions. On such a basis, for the first time oscillatory 
dependence of the penetrability on the localization of start of the packet and resonant levels of the energy of 
radiation E(ja,d) (where the penetration extremely increases) are opened (which are hidden in the semiclassical 
picture). In finishing, in Sec. 0] conclusions of the obtained results are formulated. 



2 Theoretical approach 

2.1 A model in the Priedmann— Robertson— Walker metric with radiation and gen- 
eralized Chaplygin gas 

We shall begin from consideration of a closed (fc = 1) FRW model in presence of a positive cosmological 
constant A > 0, radiation and the component of Chaplygin gas. Let us choose a minisuperspace Lagrangian in 
the following form (see Ref. [24] ) : 

where a is scale factor, a is derivative of a with respect to time coordinate t, p (a) is a general expression for the 
energy density, prad (a) is component describing the radiation in the initial stage (equation of state for radiation 
is p (a) = Prad(a)/3, p is pressure), a is parameter of Chaplygin gas (for details, see Refs. [5BH1S], also historical 
paper [29] ) . The passage to the quantum description of the evolution of the Universe is obtained by the standard 
procedure of canonical quantization in the Dirac formalism for systems with constraints. In result, we obtain 
the Wheeler~De Witt (WDW) equation (see Ref. 8^, also [TMT5] ). which after multiplication on factor and 
passage of the item with radiation p^ad to right part transforms into the following form (see Ref. .24!): 



Tr^ + V{a)\ ip{a) = E.^^d <^(a), -Brad = tt-^ 
oa'^ I 27r G 



^ "^12 + ^ ('*) r '^(") ^ ^i-ad (^(a)' -Brad = 7^777- (2) 



where 

y(a)=(-i-) ka^--^aHpj, + ^^] (3) 

and f{a) is wave function of Universe. This equation looks similar to one-dimensional stationary Schrodinger 
equation on semiaxis (of variable a) at energy i^i-ad with potential V (a) . It is convenient to use system of units 
where SttG = Af~^ = 1, and to rewrite V (a) in a generalized form as 

/ ni t \ 1/(1+") 

Via) = Aa^-Ba^{A+J^) , (4) 

where A and B are constants. In particular, at large a and A — 36, i? = 12 A this potential coincides with fl9|. 
In order to estimate ability of the approach developed in this paper below, for comparative analysis let us 
use results in [20] where a non-stationary case of the WDW equation 

^ ^ - K=ff (a)) ^{a, r)^^^^ *(a, r) (5) 

with the potential for the closed FRW model with the included generalized Chaplygin gas 

V,^{a)^3a'-^^A+^ (6) 

was studied. After change of variable a^cw = «oid V^ stationary limit of eq. ([5]) transforms into our eq. ([2]) as 
the Vef[ potential is independent on the r variable (such a choice keeps correspondence between energy levels 
that is convenient in comparative analysis). The potential (jS]) after such a transformation is shown in Fig. [1] 
and we shall study behavior of the wave function concerning it. 

2.2 Tunnehng of the packet through the rectangular barrier 

Before further solution of our problem, let us consider a general problem of quantum tunneling of a packet 
through the barrier used in cosmological models. We shall study such a process on the basis of developed 
formalism of multiple internal reflections (see main ideas, formalism, proof for tunneling, analysis of peculiarities 
of such approach in Refs. [?, ?, 25,30-33 , see also the first papers (36H38] . some development and applications 
in Refs. [39]). But, at first, we shall start from consideration of main idea of the multiple internal refiections in 
description of tunneling of the packet on the positive semiaxis of the scale factor a where possible oscillations of 
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Figure 1: Behavior of tiie potential (|6]) after ciiange anew = floid \/l2 at ^ = 0.001 and B = 0.001: (a) siiape of tiie 
barrier (Kiax = 223.52 at a = 42.322); (b) tliere is a little internal well close to zero (Knin = —8.44 at a = 0.00581) 



such a packet should be included inside internal well and, so, we shall choose the simplest potential V{a) [24:: 
V{a) = — Vb for < a < i?i (internal region I), V{a) — V\ for Ri < a < R2 (region II of the barrier) and 
V{a) = for a > i?2 (external region III). For simplicity, let us analyze a case when total energy of system E 
is higher then the barrier height Vi: E > Vi. 

In the first step we consider start of the packet in the region I which is assumed to propagate to the right 
and is incident on the first boundary of the barrier at Ri : 



^W(a,i) = / g{E - E) e''"'"^'^'/'' dE at < a < i?i. 



(7) 



where fci = y^E + Vb, E is the energy. The weight amplitude g{E — E) can be used in standard form of gaussian 
and satisfies to normalization / \g{E — E)\'^ dE = 1, value E is an average energy. This packet transforms into 
two new packets: the first packet transmitted through this boundary and propagating further in the region II, 
and the second one reflected from the boundary and propagating back in the region I: 






J g{E- E)a^^'> e""^"-'^*/^ dE at i?i < a < i?2. 



(8) 



/ g{E - E) A-^' e-^^^''-^^*'^- dE at < a < i?i, 



where ^2 = \/E — Vi . We find new unknown coefficients a^^' and Aj^ , using requirements of continuity of the 
total wave function ip(a,t) (which is summation of all packets) and its derivative at Rr. 



v(l) = 



2fci 

fci + fc2 



j,i (fci— fc2) -Ri 



(1) _ fci -h ^2»fcifli 



fci 



(9) 
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In the second step we consider further propagation of the packet iplj. {a,t), which is incident on the second 
boundary at R2. It transforms into two new packets: the first packet transmitted through this boundary and 
propagating in the region III, and the second one reflected from this boundary and propagating back in the 
region II. We define these packets in the form 



V'(2)(a,i) 



/ g{E^E)'f(^^a)e-'^*/^ dE, 



where as the stationary parts we use: 

'/'inc(a) = a^^^e**^^'^ at Ri < a < R2, 
''^'^"''"' ata>i?2, 

at Ri < a < R2, 



<'ia) = A. 



(10) 



(11) 



-ik2(i 



where k = \/ E. Imposing condition of continuity on the total wave function and its derivative at i?2, we obtain 
two new equations, from which we find new unknowns coefficients A>p and ^^^': 



/3(i)=i?+.a(i), R\ 



kV-l ^. ^ (12) 



•■2 ■ " ' ^ '■2 ~ • • ^ 



2i fe2 fe 



k2 + k 

We have introduced two new coefficients T^ , R^ , which logically connect the transmitted and reflected ampli- 
tudes Arp and /3*^^-' with the incident amplitude a*-^^ in this step. Here we shall use bottom index for denotation 
of number of the considered boundary, upper (top) sign "+" or "— " for positive (to the right) or negative (to 

the left) direction of the incident wave, correspondingly. So, we can write T^ = a^^'' and Ri — A)^' also. 

(2) 
In the third step we consider further propagation of the reflected packet tp^J in the region II. Incidenting 

(3) 
on the first boundary, it transforms into new packet ■0jj. , transmitted through this boundary and propagating 

in the region I, and into new packet i/'rof ' reflected from boundary and propagating back in the region II. We 
define the new packets by eq. (jlOp (with upper index 3), where as the stationary parts we use: 

(13) 



^l;^l(«) ^ 


= ^L1(«) 


a.t Ri < a < R2 


Af(a) - 


_ 4(2) -ifeia 


at < a < i?i , 


^i^^(«) = 


= Q,(2)e«fe2a 


at i?i < a < i?2 



From continuity conditions for the total wave function and its derivative at i?i we find the unknowns coefficients 
Ag^ anda(2). 

^2 



' ^' '^ -' lltll ^.. „ (14) 



^ ^ ki+k2 

In the forth step we need to consider further propagation of the reflected packet ^pj.J in the region I in 
the 1-st step. It is incident on the first boundary at a = transforming into new packet propagated to the 
right. At such a point we can include different considerations of origin of possible sources at a = 0, possible full 
propagation (like in spherically symmetric problems of quantum decay in nuclear physics which is 3-dimensional 
and we have no additional boundaries at a = 0) or, in contrary, possible full refiection used in different quantum 
approaches (like introduction of an infinite potential wall at a = in Ref. [19 ). In order to produce ability to 
work with different such considerations, we write: 

^!nc(«)-^rcf(«), At\a,k^) = R^-^[%,^k,)=4tU^''^'^ atO<a<i?i, (15) 

where 

Supposing full propagation through this boundary (without any possible reflections), we obtain Rq = — 1. If 
we liked to include the infinite potential wall at a = 0, than we should have Rq = —1 also. 

Analyzing further reflections and transmission of the packets concerning the boundaries by such a way, 
we have concluded that any following step is similar to one from 4 considered above. From analysis of these 
steps we have found recurrent relations for calculation of unknown amplitudes ^[nj, A^^, AIj!'' a^"^ and /S^"^ 
for arbitrary step n, and we have calculated summations of these amplitudes. However, these series could be 
calculated easier, if to apply coefficients T^ and Rf. Analyzing all possible "paths" of the propagations of all 
possible packets inside the barrier and internal well, we obtain: 

E A^l = i + RtRo+RtRo-RtRo+- = ^+ E {Rt Ro ' 



A^= 1^-42- (16) 



n=l 

-t-00 , ^ /+00 , X \ /- -1 /+OC 



1 — Ri Rq 



'2 ' 



(17) 



«f . (1 + E («.- «f )■•■) - -^ - (E ^L1) - «f 



I- Rq Ri 



where 

Rt = Rt + r+ Rt T,- . (i + E {Rt Ri)"') =Rt+ f f/'l , 

, +00 -\ T^+ '-'''''' (18) 



m=l ' 1 — i?2 ^1 

The resultant expressions for the incident, transmitted and reflected packets concerning the barrier are 
written in form of eq. ()10p . where the following stationary wave functions should be used: 



¥'inc(a) = f^ikia^ forO<a<i?i, 

+00 
(ptr(a) = E ^re''^"' ^r a > i?2, 

+ CX3 

'y'rofla) = E ^fle-*'^!'', forO<a<i?i. 

n=0 



(19) 



At finishing, we determine the full amplitudes 

AT-E4"^ An^'^Af, «=E"("'-S, /3= E/3'"^=«-^2+ (20) 

n— 1 n— 1 n— 1 -^2 n— 1 

and coefficients T and i? describing penetration of the packet from the internal region outside and its reflection 
from the barrier 

TmIR= -r~ \At\ — L4inc • Jbar, Tbar = 1~ K 1 i , , 

Kl ' ' ' ' fci ' ' (21) 

RmIR = \A-r\ — |Ainc| • ^bar, -Rbar = |i?i | , 

where Tbar and i?bar are coefficients of penetrability and reflection of the barrier (in standard definition), and 
I 1 2 
Amc is coefficient determining oscillations of the packet inside the internal region (this is fully quantum analog 

of the normalization factor F introduced in Ref. |40j for semiclassical description of nuclear decay) . 

Series E^inc S ^t ' TIi^r ^ E o'^"'' ^"^^ S /^^"^ obtained using the approach of the multiple internal 
reflections, exactly coincide with the corresponding coefficients Ainc, At, Ar, a and (i calculated by standard 
stationary method (where the continuity conditions of the stationary total wave function and its derivative are 
used at each boundaries, and the wave function is not equal to zero at a — Q). We test property: 

— \At\'' + \Ar\^ = 1 or Tmir + Rmir^I, (22) 

ki 

which is fulfilled and confirms that the method MIR gives us proper solution for the wave function. If energy is 
less then the height of the barrier, then for description of penetration of the wave through such a barrier with 
its tunneling it needs to use the following change [25 11311133] : 



fe^i?, ^=x/F"^TT. (23) 

Using it, all found above solutions are applied for the problem with tunneling through the barrier. 

2.3 Tunneling of the packet through barrier composed from arbitrary number of 
rectangular steps 

Now let us come to another essentially more difficult problem of the packet penetrating through the radial 
barrier of arbitrary shape in the cosmological problem. In order to apply the idea of multiple internal refections 
for study the packet tunneling through the real barrier, we have to generalize formalism of the multiple internal 
reflections presented above. We shall assume that the total potential has successfully been approximated by 
finite number N of rectangular steps: 

{Vi, at flinin < a < fli (region 1), 

V2, at ai < a < a2 (region 2), .^^. 

Vn, at QN-i < a< flniax (region N), 

where Vi are constants (i = 1 . . . N) . Now let us assume that the packet starts to propagate outside inside the 
region with some arbitrary number M (for simplicity, we denote its left boundary ajvf-i as astart) from the left 



if [a) 



(25) 



of the barrier. We shall be interesting in solutions for above barrier energies while the solution for tunneling 
could be obtained after by change i ^i — > fc^ . A general solution of the wave function (up to its normalization) 
has the following form: 

ai e"-^^"- + Pi e-*'=i'^, at amin < a < ai (region 1), 

ftM-i e*'^"-!'^ + pM-i e-''^■*'-l^ at aM-2 <a< au-i (region M - 1), 
^ikMa _^ ^^ g-jfcMa^ ^1- ^^_^ < a < flA/ (region M), 

QAf+i e*'=^^'+i'^ + /3m+i e-'^^M+^a^ ^^ aM <a< oaz+i (region M + 1), 

a„_ie*''"-i" + /3Ar_ie-''=«-i", at aAr_2 < a < ajv_i (region TV - 1), 
At e*''"", at oat-i < a < a,nax (region N), 

where aj and /Sj are unknown amplitudes, At and An are unknown amplitudes of transmission and reflection, 
ki = j^^j2m{E — Vi) are complex wave numbers. We have fixed the normalization so that modulus of the 
starting wave e**^"" equals to one. We shall be looking for solution for such a problem by the approach of the 
multiple internal reflections. 

Let us begin from consideration of start of the packet in the region with number M , which propagates to 
the right. At first, we shall study its propagation inside the right part of the potential with barrier, starting 
from this region. According to the method of the multiple internal reflections, scattering of the packet on 
the barrier is considered consequently by steps of its propagation relatively to each boundary of the barrier 
(the most clearly idea of such approach can be understood in the problem of tunneling through the simplest 
rectangular barrier, see [25l[3TJ[33] where one can find proof of this fully quantum exactly solvable method, one 
can analyze its properties) . Each step in such consideration of propagation of the packet will be similar to one 
from the first 2N — 1 steps, independent between themselves. From analysis of these steps recurrent relations 
are found for calculation of unknown amplitudes A^"'^ S'("\ a^") and /3^"^ for arbitrary step n, summation of 
these amplitudes are calculated. We shall be looking for the unknown amplitudes, requiring wave function and 
its derivative to be continuous at each boundary. We shall consider the coefficients T^ , T^ ■ ■ ■ and Rf , R^ 
. . . as additional factors to amplitudes e***^". Here, bottom index denotes number of the region, upper (top) 
signs "+" and "— " denote directions of the wave to the right or to the left, correspondingly. At the first, we 
calculate T^ , T^ ■ ■ ■ T^_i and Rf , K^ ... R^-i- 

Kj + Kj-f-i Kj + fcj + 1 , , 

(26) 

, _ kj — kj+l 2ikjaj p- _ fcj + 1 — kj 2ikj + iaj 



] 



j ' :y+i j ' i+i 



Analyzing all possible "paths" of the propagations of all possible packets inside the barrier and internal well, 
we obtain: 



n— 1 m—1 1 RaI M — 1 

+ 00 , . /+00 



E ^T ~ ( E ^inc j ' 1 -'w-2 -'- N-1 + -'^ N-2 ' ^N-1 ^N-2 ' -'- N-1 + ■•• f ~ I E ^inc j ' ^J 

1=1 ^n=l ' '- -' ^n=l ^ 

fcx) ,, ^ ^_ ^ ^^ __ 

E "^R. — ^M + ^M ' ^M-l ^M + ^M ' ^M-1 ^M ' ^M-1 ^M + ••■ = 

1=1 

= -Rm • (^1 + E {Rm-1 Rai)"') = __ ^J' n+ = [z2 Anc) • Rtl 



N-1' 



(27) 



where 



^j-i - 



^ rn=l ^ 1 - RJRj_i 



>-M' 



TU^T^i 



+ OC' 



^7+1 = ^7+1 + T^iR^T^i (1 + E (^;+i^7)"0 = ^7+1 + ^^""'5 It- - (28) 



^i = ^^i(i+ E (i?^iiV^'"^ - ^^' 



1 — Rj_^iR 



j+l-"-] 



and selecting as starting the following values: 

we calculate successively coefficients R%_2 ■ ■ ■ RXi^ ^m+i ■ ■ ■ ^n-i ^^"^ ^aj+i • • ■ '^n-v 

At second, we shall consider further propagation of all packets which propagate in the region with number M 
to the left. Such packets are formed in result of all possible reflections from the right part of potential, starting 
from the boundary oa/. In the previous section for description of their reflection from the left boundary i?o to 
the right we used coefficient R^ . Now if we liked to pass from simple boundary aM~i to the left part of the 
potential well starting from this point up to Omin, we should generalize the coefficient i?^/-! '-"^ ^m-i- -"-^ turns 
out that the middle recurrent formula p8|) is absolutely applicable in such a case also, where for definition of 



(30) 



Tj and R^ we should use eqs. ([26|) again. At finishing, we determine coefficients aj and 13 j: 

— „(") — T+ /I I Y^ /'P+P- \m\ _ J-1 _ _J_ 

- R^RJ^i T+ ' 



1=1 ^ m=l ^ 1 



E pf = TU 1 + E iRtRj-^r) Rt - , %,"- = -^ 

n=l ^ m=l ' 1 - RJR,_. T^ 



"-J '-''j-1 3 



the amplitudes of transmission and reflection: 



+ 00 , +00 , , +00 , 2^+ +00 , 



n—l n—1 n—1 " i- j n—1 



and coefficients T and R describing penetration of the packet from the internal region outside and its reflection 
from the barrier 

T^ - ^^ U |2 I /I |2 T^ T ^N \^+ \2 
^MIR=- \At\ = |Ainc| -Jbar, J bar = 1 I^AT-l] J 



A/ ' ' ' ' ku I '^ ^1 (32) 

+ r 

Af I 



RaIIR — \^B.\ — \M-ac\ ■ Rban Rbm — \R 



Choosing amin = 0, we assume full propagation of the packet through such a boundary (without any possible 
reflection) and we have Rq = — 1 (it could be interesting to analyze results varying Rq). We check the property: 

^\At\^ + \Ar\^ = 1 or Tmir + Rmib.^1, (33) 

Km 

which should be the test, whether the method MIR gives us proper solution for wave function. Now if energy 
of the packet is located below then height of one step with number m, then for description of transition of this 
packet through such barrier with its tunneling it shall need to use the following change: 

fcm->iCm- (34) 

For the barrier consisting from two rectangular steps of arbitrary heights and widths we have already obtained 
coincidence between amplitudes calculated by method of MIR and the corresponding amplitudes found by 
standard approach of quantum mechanics up to first 15 digits. Increasing of number of steps up to some 
thousands keeps such accuracy and fulfillment of the property ([33]) (see Appendix |X] where we present shortly 
the standard technique of quantum mechanics applied for the potential p4| and all obtained amplitudes). This 
is important test which confirms reliability of the method MIR. So, we have obtained full coincidence between 
all amplitudes, calculated by method MIR and by standard approach of quantum mechanics, and that is way 
we generalize the method MIR for description of tunneling of the packet through potential, consisting from 
arbitrary number of rectangular barriers and wells of arbitrary sizes. 

3 Results 

We have applied the method above to analysis of behavior of the packet in its propagation relatively the barrier 
(|6|) (taking into account Onow —> vl2 Ooid)- The first interesting results which we have obtained is visible change 
of the penetrability on displacement of the starting point Raiin ^ f l£ Ri, where we putted the packet for start. 
Using possibility to decrease width of intervals up to enough small limit (and choosing, for conveniens, the 



width of each interval to be the same) , we call i?miii starting point, from where the packet begins its propagation 
outside. We have analyzed how much position of such a point influences on the penetrability. In Fig. [5] one 
can see that at arbitrary fixed energy of radiation -Erad the penetrability of the barrier is changed strongly in 
dependence on -Rmin: it has oscillating behavior, difference between its minimums and maximums is minimal 
at i?min in the center of the well (i. e. its change tends to zero in the center of the well), with increasing of 
Rmin this difference increases, achieving to the maximum close to the turning point. On such a basis we (for 
the first time) establish dependence of penetrability on the starting point Rstnit of the packet. Behavior of the 
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Figure 2: Dependencies of the coefficients of penetrability, reflection, oscillations and penetration on the position of the 
starting point, Oatart for the energy E — 220 {A = 0.001, B — 0.001, amax = 70, total number of intervals is 2000): (a) 
penetrability; (b) reflection; (c) coefficient of oscillation; (d) coefficient of penetration. For all presented values we have 
achieved accuracy [Tbar + Rhm — 1] < 1 ■ 10~^^. 



coefficients of reflection, oscillations and penetration turn out to be similar. 

Usually, in cosmological quantum models the penetrability is determined by the barrier shape. While in 
non-stationary approach one can flnd papers where a role of the initial condition is analyzed in calculations of 
rates, penetrability etc. (such papers are very rare and questions about dynamics have not been studied deeply), 
then stationary limit does not give us a choice for work, practically. Note that such an understanding about the 
penetrability of the barrier is prevailing till present day not only in quantum cosmology, but also in problems 
of decays and fission in theoretical nuclear physics. In last topic, such a point of view is as much deep-rooted 
as description of dynamics of decays and fission, calculations of half-lives are fulfilled without inclusion of such 
initial condition. Let us give only some examples. In Ref. 41] agreement between experimental data of a-decay 
half- lives and ones calculated by theory is demonstrated in a wide region of nuclei from ^°^Te up to nucleus 
with Ad = 266 and Zd = 109 (see Ref. "^^jfiS] for improved approaches and data of last years). In review [33] 
methodology of calculation of half-lives for spontaneous-fission is presented (see eqs. (21)-(24) in p. 321). In 
prevailing two-potential approach (TPA) and semiclassical approach (WKBA) for determination of half-lives of 
nuclear decay by emission of proton (nuclear proton decay), which today are recognized as the most accurate 



and reliable (for example, see Refs. [iniH51 - H7] ) . influence of such initial condition on results has not been taken 
into account and has not been analyzed yet (for details and explanations, see [33]). Such approaches forms 
grounds for producing tables of nuclear data (for example, see [H])- The unexpected result in Fig. [5] leads to 
inevitable change of understanding about the penetrability. The first and direct conclusion is: the penetrability 
should be connected with the initial condition (not only in non-stationary consideration, but also in stationary 
one), which sets localization of start of the packet outside. The second important conclusion is: even in the 
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Figure 3: Dependencies of the coefficients of penetrability, reflection, oscillations and penetration on the position of the 
external region, flmax for the energy E = 223 {A = 0.001, B = 0.001): (a) penetrability; (b) reflection; (c) coefiicient of 
oscillation; (d) coefficient of penetration. For all presented values we have achieved accuracy |Tbar + Rha.T — 1| < 1 ■ 10~^^ 
(majdmum number of intervals is 2000). 

stationary consideration the penetrability of the barrier should be determined in dependence on the initial 
condition. 

The first natural question, which could appear from analysis of such results, is how much they are reliable. In 
particular, would such results be destroyed if we shifted the external boundary outside (while in the semiclassical 
case we are restricted by two turning points only, then in the fully quantum approach the external tail of the 
barrier effects on the results inevitably and additionally)? Taking strong decreasing of the external tail of the 
barrier into account to minus infinity, one could even expect for this. Results of such calculations are presented 
in Fig. [21 where it is shown how the penetrability is changed at increasing the external boundary Omax (for the 
most clearness, we have fixed the starting point: Ogtart = 10)- One can see that all calculations are perfectly 
convergent, that confirms efficiency of the techniques of the multiple internal reflections. This points to that 
difficult attempts and aspiration to determine the outgoing wave in the asymptotic limit properly (on which the 
boundary condition of tunneling is based) have no any practical sense. Therefore, we have chosen Omax = 70 
for further calculations with the studied potential. However, one can see that inclusion of the external region 
in calculations can change the coefficients of penetrability and penetration up to 2 times for the chosen energy 
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Figure 4: Inside the energy region -Erad = 200 — 223 we have observed 19 resonant levels (we have choose: astart = 10, 
flmax ~ 70, number of intervals inside the scale axis a is 1000, number of intervals of energy is 100000): (a) the coefflcient 
of penetration in dependence of the -Erad energy; (b) the penetrability in dependence of the Erad energy; (c) the coefficient 
of oscillation in dependence of the E^rad energy; (d) difference between two closest resonances Eres, next — Erea, previous in 
dependence on the Eros energy. 

level. In Tabl. [1] one can see how the coefficients of the penetrability and reflection are changed for different 
values of the starting point for the same energy. One can see that, in contrast to other known approaches, our 
developed method in determination of these coefficients gives the fantastic accuracy: |Tbar + ^bar — 1| < 1-10~^^. 

The second question of not less importance is how large this effect in calculations of the penetrability: if 
it was small than, for example, the semiclassical approaches would have enough good approximation. From 
Tabl. [T] it follows that the penetrability is changed not strongly in dependence on shift of the starting point. 
However, such small variations are connected with relatively small height of the barrier and depth of the well, 
while they would be not small at another choice of parameters (the coefficient of oscillation and penetration 
turn out to be changed essentially at some definite energies of radiation, see below). So, this effect is supposed 
to be larger at increasing height of the barrier and depth of the well, and also for near-barrier (i. e. for energies 
comparable with the barrier height) and above-barrier energies of radiation. 

Then, we have analyzed how these characteristics are changed in dependence on the energy of radiation. 
Result has turned out to be over- unexpected (see Fig. |4]): the coefficient of penetration has oscillating behavior 
where peaks are clearly shown, located about the same distances, between which smooth and stable in calcula- 
tions wells are observed with minimums! By other words, in the fully quantum approach we for the first time 
have observed clear pictures of resonances, which could be connected with quasi-stationary states. More detailed 
analysis gives: at increasing energy of radiation the penetrability is changed monotonously and determines a 
general tendency of change of the coefficient of penetration, while the coefficient of oscillations introduces the 
peaks. Now a reason of presence of resonances has become clearer: oscillations of the packet inside the internal 
well cause them, while the possibility of the packet to penetrate through the barrier (described by the penc- 
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trability of the barrier) does not take influence on them absolutely. In Tabl. [2] we present the resonant levels 
calculated by such a way inside the energy range E'rad = 200-223, and the coefficients connected with them. 
To such data we add a full list of the found resonant levels inside energy range i^i-ad = 0-200 in Tabl. [3] (we 
have established 134 resonant levels). We seem for the first time to have obtained the quasi-stationary levels in 
the problem of decay in the fully quantum approach, separating the coefficient of penetration (having sense of 
width in theory of quantum scattering and decays) on the coefficients of penetrability and oscillations. 

4 Conclusions 

The new method for determination of probability of penetration of the packet from the internal well outside 
with its tunneling through one-dimensional barrier of arbitrary shape used in problems of quantum cosmology, 
is presented. Note the following: 

1. The method is further development of approach of multiple internal reflections (see Refs. [?. ?.[25l[30l- 
[55] . also Refs. [55H5^ ). where a process of tunneling of the packet through the barrier is considered 
consequently by steps of its propagation relatively to each boundary of the barrier. The method is fully 
quantum, allows to determine amplitudes of wave function, coefficients of penetrability and reflection 
relatively the barrier. For the first time exact analytical solutions for amplitudes of the wave function, 
penetrability T and reflection R for the barrier, composed from arbitrary number n of rectangular potential 
steps, are found. At limit n — ;> oo these solutions could be considered as exact limits for potential with 
interested barrier and internal well of arbitrary shape. 

2. Accuracy of the method in determination of penetrability Tbar and reflection /?bar is: iTbar + ^bar ~ 1| < 
1 • 10~^^ (see Tabl.[T]). Author has not found other methods achieving such accuracy in similar problems of 
quantum physics (with possible exception of some selected cases of exactly solvable barriers which could 
be obtained by methods of supersymmery). 

3. On the basis of the method the probability of penetration of the packet from the internal well outside 
with its tunneling through the barrier of arbitrary shape called coefficient of penetration, is determined. It 
succeeds to separate that coefficient explicitly on the penetrability and new coefficient, which characterizes 
oscillating behavior of the packet inside the internal well and is called coefficient of oscillations. That 
found for the first time formula seems to be fully quantum analogue of the semiclassical formula of F width 
of decay in quasistationary state proposed by Gurvitz and Kalbermann in Ref. |40| (here, the coefficient 
of oscillations is fully quantum analogue for the semiclassical F factor of formation and the coefficient of 
penetration is analogue for the semiclassical F width). 

This method has been applied for study of properties of the packet, describing evolution of universe on the 
first stage in the closed Friedmann-Robertson- Walker model with quantization in the presence of the positive 
cosmological constant, radiation and component of generalize Chaplygin gas with potential chosen from [20]. 
Let us formulate main results obtained for the first time: 

1. For the same chosen energy of radiation Erad the penetrability of the barrier is changed visibly in de- 
pendence on the position of the starting point i?start inside the internal well, where the packet begins to 
propagate (see Fig. [2|): the penetrability has oscillating behavior, difference between its minimums and 
maximums is minimal at i?start in the center of the well, with increasing i?start this difference increases 
achieving maximum close to the turning point. The behavior of the coefficients of reflection, oscillations 
and penetration turns out to be similar. Coincidence (up to the first 15 digits) of the amplitudes of 
the wave function obtained by such a method, with corresponding amplitudes obtained in the standard 
approach of quantum mechanics (see App. EJ at different energies E^rad confirms that this result does not 
depend on a choice of the fully quantum method applied for calculations. Such a peculiarity is shown 
in the fully quantum non-stationary and stationary considerations and it is hidden after imposing the 
semiclassical restrictions. 

2. In non-stationary and stationary considerations the penetrability of the barrier should be connected with 
initial condition localizing start of the packet. Note that possible introduction of the initial condition into 
known stationary semiclassical models could change their results. 

3. The penetrability is changed visibly, if to take the external tail of the barrier into account. For example, 
for the barrier ([6]) with parameters A = 0.001 and B = 0.001 (see Fig. [1]) at the energy of radiation 
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-E-rad = 223 the penetrability is changed up to 2 times (see Fig. ^. If to increase the external boundary 
flmax , all amplitudes and coefficients are convergent in calculations that confirms efficiency of the developed 
method. 

4. Dependence of the coefficient of penetration on the energy of radiation has oscillating behavior: here peaks 
are clearly shown, localized approximately on the same distances, between which smooth minimums (wells) 
stable in calculations are observed (see Fig. |4]). By other words, for the first time in the fully quantum 
approach we have obtained clear and stable picture of resonances, which indicate on presence of some 
quasistationary states. Here, with increasing of the energy of radiation the penetrability is changed 
monotonously and it describes a general tendency of behavior of the coefficient of penetration, while the 
coefficient of oscillations gives peaks. Now a reason of existence of resonances becomes clear — oscillations 
of the packet inside the internal well give them. For example, for the barrier (jB)) with parameters A = 0.001 
and B = 0.001 I establish 134 such resonant levels inside range E^ad — 0-223 (see Tabl. [51 [3]) . 

5. The dependence of the penetrability on the starting point has clear established maximums and minimums. 
On such a basis one can suppose that the most probable start of the packet (describing start of expansion 
of the universe) is in one point of such maximums. This allows to predict some definite initial values of 
the scale factor, at which the universe begins to expand (such initial data is direct result of quantization 
of the classical cosmological model) . 

6. Modulus of the wave function both in the internal, and in the external regions has clear established own 
maximums and minimums |24ll49j . This indicates on such values of the scale factor, at which "appearance" 
of the universe will be more or less probable. By other words, radius of the universe during its expansion 
is changed not continuously, but passes consequently through some definite discrete values connected 
with these maximums. It follows that after quantization space-time of universe on the first stage of its 
expansion seems to be rather discrete than continuous. According to results [24t49 . difference between 
maximums and minimums with increasing of the scale factor a is slowly smoothed and we obtain obvious 
for us continuous structure of the space-time at latter times. Discontinuity of space-time is direct result 
of quantization of cosmological model, which is shown the most strongly on the ffist stage of expansion 
and disappeared after imposition of the semiclassical approximations. 

A Direct method 

We shall add shortly solution for amplitudes of the wave function obtained by standard technique of quantum 
mechanics which could be obtained if to use only condition of continuity of the wave function and its derivative 
at each boundary, but on the whole region of the studied potential. At first, we find functions /2 and 52 (from 
the first boundary): 

f^ ^ h±A e^'''-^\ g^ = J^e'i^+^-^-K (35) 

k2 — k k — k2 

Then, using the following recurrent relations: 

-• ^ (fcj + i - kj) e^'^^"^- -f- f, {kj + i + kj) 2»fe, + ix, /of.^ 

^■'+' (fc,+i+fc,)e2''=.-.+/,(fc,+i-A:,)'^ ' ^^^^ 



we calculate next functions f^, fi, /s ■ ■ ■ fn, and by such a formula: 

2k. gi(kj + l+kj)xj 

(fc,+i + /c,)e2^fe.-. +/,(fcj+i - kj) 
the functions 53, (74, g^ ■■■gn- From /„ and g^ we find amplitudes a„, /?„ and amplitude of transmission Aq^: 






/3„ = 0, Ar^an^-^. (38) 



Now using the recurrent relations: 



^ _ pikjXj--i \ P.p ikjXj — i — „. „ %kj-iXj-i 

«j-i = ^^ T^^^ 7 ^^^ (39) 
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and such a formula: 

ft— «.-•/. +5,, (40) 

we consistently calculate the amplitudes a„_i, /3„_i, a„_2, /3n-2 • • ■Ck2, /?2- At finishing, we find amplitude of 
reflection Afj/. 

An = a2 e'('=+'=2)^i + /32 e'^*^"*^^)"^! - e^*''"'^ . (4f ) 

As test we use condition: 

t^\AT\' + \An\'^l. (42) 

fci 

In order to check all amplitudes obtained previously by the MIR approach, we have used such a techniques 
and obtained coincidence up to the first 15 digits for all considered amplitudes. In particular, we reconstruct 
completely the pictures of the probability and reflection presented in Fig. [5] (a) and (b). Fig. |3] (a) and (b). 
Fig. 131(b), but using standard technique above. So, result on the oscillating dependence of the penetrability of 
the position of the starting point Riorm i*^ such figures is independent on the fully quantum method chosen for 
calculations. 
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0.289642192688668 
0.289696448926345 
0.289644249841587 
0.289637485908461 
0.289672548079996 
0.289672862933907 
0.289611311229824 
0.289731953026223 
0.289575702298831 
0.289737152195946 
0.289552884375544 
0.289669677006127 
0.289791629591448 
0.289578213304735 
0.289545256222507 
0.289581718367642 
0.289745254914410 
0.289708963481124 
0.289846866037466 
0.289824379615929 
0.289408078297668 
0.289904994878596 
0.289330571976127 
0.289506924291343 
0.289723176821301 
0.290117973838998 
0.289691943875087 
0.289084110255789 
0.290388261985181 
0.288817117691736 
0.289619296892063 
0.290071707190809 
0.289193475557903 
0.287841171185186 
0.291902880497247 
0.286694489516019 
0.287618039346707 
0.285142815727187 
0.279635849388475 
0.329388280906961 
0.292960219409352 
0.353172326894230 
0.641964142658414 



0.0252780829608057 
0.0252782910350983 
0.0252776048415682 
0.0252773925052633 
0.0252771781182692 
0.0252769611361457 
0.0252794996594910 
0.0252782017545698 
0.0252753096406404 
0.0252805857676684 
0.0252792667386222 
0.0252806184777316 
0.0252811365836465 
0.0252723581940758 
0.0252849928275708 
0.0252707897908712 
0.0252693323305645 
0.0252868210542670 
0.0252726558228249 
0.0252731048676081 
0.0252714868761523 
0.0252620108618688 
0.0252784521274925 
0.0252899361825347 
0.0252871585222866 
0.0252953619874030 
0.0253069628091520 
0.0253096065275343 
0.0252656032511840 
0.0253120503626760 
0.0252771291479563 
0.0252663105653185 
0.0252687898567804 
0.0252237178217930 
0.0252252761004983 
0.0252787062232434 
0.0252066359105951 
0.0252789797198770 
0.0253000826894572 
0.0252162385846372 
0.0251410123308979 
0.0251470261390868 
0.0250645870558445 
0.0251150672964445 
0.0255534898499390 
0.0252959368485190 
0.0245869086911412 
0.0242534272927471 
0.0256960991603236 
0.0267084356304886 
0.0443847180924309 
0.0501465994586273 
0.0694934432562353 



Table 1: Dependencies of the coefficients of penetrabihty Tbar, reflection -Rbar, oscillations A'osc and penetration Tmir. 
on the position of the starting point flstart for the energy E — 223 (we used A — 0.001, B — 0.001). For all presented 
values we have achieved accuracy |Tbar + Rhur — Ij < 1 ■ 10^^''. 
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Energy, 


Fully quantum method 


A Energy, 

-C^rad,ncxt -C/j-ad, prev 


Penetrability, Tbar 


Oscillation, Kqsc 


Penetration, T^im 


200.631586315863 
201.993199931999 
203.342393423934 
204.678706787068 
206.001450014500 
207.310163101631 
208.603926039260 
209.882278822788 
211.144301443014 
212.389073890739 
213.614986149861 
214.821118211182 
216.005400054000 
217.166221662217 
218.300363003630 
219.404604046040 
220.473424734247 
221.499004990050 
222.464784647846 


2.40152454 x lO"'^^ 
1.14426158 X 10-'*2 
5.14289866 X 10"*° 
2.15097704 X lO'^^ 
8.44746902 x lO'^^ 
3.08154101 X 10-32 
1.04396547 x lO^^g 
3.28814440 x lO'^^ 
9.50957150 x lO^^s 
2.54399190 x 10^22 
6.19406186 X lO^^o 
1.37175187 X 10-1^ 
2.75087788 x lO"*'^ 
4.91371044 X lO^^^ 
7.77746865 x 10"" 
1.07339618 X 10-^ 
1.25345583 x 10"^ 
1.20673982 x 10"* 
8.80196772 x 10"^ 


1.64916505 X 10+^ 
1.89183111 X 10+*^ 
1.45030359 x 10+«^ 
6.62174746 x 10+^ 
6.74278273 x 10+^ 
4.64608634 x 10+^ 
1.62095089 x 10+*^ 
2.49850310 x 10+^ 
1.65795995 x 10+^ 
3.64786127 x 10+^ 
9.03599338 x 10+^ 
5.19487954 x 10+^ 
2.15778829 x 10+^ 
4.16964758 x 10+^ 
7.48378112 x 10+^ 
3.58886110 X 10+^ 
8.94224185 x 10+^ 
3.66509631 x 10+^ 
4.20824526 x 10+* 


3.96051035 x lO'^^ 
2.16474967 X lO'^^ 
7.45876441 x 10-^4 
1.42432267 x lO'^i 
5.69594482 x lO'^^ 
1.43171056 X 10-26 
1.69221676 x lO^^s 
8.21543901 X lO'^i 
1.57664887 x lO"*** 
9.28012954 x lO"*^ 
5.59695021 x 10"" 
7.12608576 X lO^^^ 
5.93581210 X 10"^ 
2.04884408 x 10"^ 
5.82048731 x 10"^ 
3.85226983 x 10"^ 
1.12087052 X 10+" 
4.42281767 X 10+* 
3.70408389 x 10+^ 


1.36161361613616 
1.34919349193491 
1.33631336313363 
1.32274322743227 
1.30871308713087 
1.29376293762937 
1.27835278352783 
1.26202262022620 
1.24477244772447 
1.22591225912259 
1.20613206132061 
1.18428184281842 
1.16082160821608 
1.13414134141341 
1.10424104241042 
1.06882068820688 
1.02558025580255 
0.96577965779657 



Table 2: The resonant levels inside the energy region -Brad = 200 - 223 (we have choose: ttstart = 10, amax = 70, 
number of intervals inside the scale axis a is 1000, number of intervals of energy is 100000). We also add the coefficients 
of penetrability Tbar, oscillations Kosc and penetration Tmir for such levels (we used A — 0.001, B — 0.001). For all 
presented energy levels, -Erad, we have achieved accuracy in calculations of Tbar and -Rbar^ |Tbar + -Rbar — 1| < 1 • 10~^^. 
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No 


Energy 


Penetrability, 


AE 




£^rad 


Tb 


ar 




1 


0.1221 


1.696 X 


10-bi8 


- 


2 


1.5261 


9.595 X 


10-512 


1.526 


3 


3.5352 


3.657 X 


10-504 


2.009 


4 


5.5377 


3.908 X 


10-497 


2.002 


5 


7.5357 


1.938 X 


10-490 


1.998 


6 


9.5293 


5.550 X 


10-484 


1.993 


7 


11.5184 


1.016 X 


10-477 


1.989 


8 


13.5031 


1.290 X 


10-471 


1.984 


9 


15.4856 


1.222 X 


10-465 


1.982 


10 


17.4614 


8.601 X 


10-460 


1.975 


11 


19.4350 


4.822 X 


10-454 


1.973 


12 


21.4042 


2.174 X 


10-448 


1.969 


13 


23.3689 


7.949 X 


10-443 


1.964 


14 


25.3314 


2.462 X 


10-437 


1.962 


15 


27.2895 


6.413 X 


10-432 


1.958 


16 


29.2431 


1.414 X 


10-426 


1.953 


17 


31.1945 


2.740 X 


10-421 


1.951 


18 


33.1414 


4.579 X 


10-416 


1.946 


19 


35.0840 


6.686 X 


10-411 


1.942 


20 


37.0243 


8.752 X 


10-406 


1.940 


21 


38.9579 


9.932 X 


10-401 


1.933 


22 


40.8893 


1.022 X 


10-395 


1.931 


23 


42.8185 


9.533 X 


10-391 


1.929 


24 


44.7410 


7.845 X 


10-386 


1.922 


25 


46.6614 


5.945 X 


10-381 


1.920 


26 


48.5772 


4.057 X 


10-376 


1.915 


27 


50.4909 


2.551 X 


10-371 


1.913 


28 


52.3979 


1.450 X 


10-366 


1.906 


29 


54.3027 


7.581 X 


10-362 


1.904 


30 


56.2030 


3.658 X 


10-357 


1.900 


31 


58.0989 


1.617 X 


10-352 


1.895 


32 


59.9926 


6.662 X 


10-348 


1.893 


33 


61.8796 


2.518 X 


10-343 


1.887 


34 


63.7644 


8.852 X 


10-339 


1.884 


35 


65.6448 


2.902 X 


10-334 


1.880 


36 


67.5207 


8.828 X 


10-330 


1.875 


37 


69.3944 


2.522 X 


10-325 


1.873 


38 


71.2614 


6.679 X 


10-321 


1.867 


39 


73.1240 


1.639 X 


10-316 


1.862 


40 


74.9844 


3.831 X 


10-312 


1.860 



No 


Energy 


Penetrability, 


AE 




Erad 


n 


ar 




41 


76.8403 


8.368 X 


10-308 


1.855 


42 


78.6918 


1.709 X 


10-303 


1.851 


43 


80.5389 


3.299 X 


10-299 


1.847 


44 


82.3793 


5.871 X 


10-295 


1.840 


45 


84.2175 


1.001 X 


10-290 


1.838 


46 


86.0534 


1.623 X 


10-286 


1.835 


47 


87.8827 


2.452 X 


10-282 


1.829 


48 


89.7076 


3.511 X 


10-278 


1.824 


49 


91.5280 


4.721 X 


10-274 


1.820 


50 


93.3440 


6.050 X 


10-270 


1.815 


51 


95.1555 


7.296 X 


10-266 


1.811 


52 


96.9626 


8.368 X 


10-262 


1.807 


53 


98.7653 


9.104 X 


10-258 


1.802 


54 


100.5613 


9.263 X 


10-254 


1.795 


55 


102.3551 


9.112 X 


10-250 


1.793 


56 


104.1444 


8.460 X 


10-246 


1.789 


57 


105.9271 


7.440 X 


10-242 


1.782 


58 


107.7075 


6.267 X 


10-238 


1.780 


59 


109.4813 


4.987 X 


10-234 


1.773 


60 


111.2507 


3.784 X 


10-230 


1.769 


61 


113.0156 


2.730 X 


10-226 


1.764 


62 


114.7739 


1.870 X 


10-222 


1.758 


63 


116.5277 


1.214 X 


10-218 


1.753 


64 


118.2771 


7.598 X 


10-215 


1.749 


65 


120.0220 


4.512 X 


10-211 


1.744 


66 


121.7625 


2.578 X 


10-207 


1.740 


67 


123.4963 


1.389 X 


10-203 


1.733 


68 


125.2257 


7.181 X 


10-200 


1.729 


69 


126.9485 


3.520 X 


10-196 


1.722 


70 


128.6668 


1.649 X 


10-192 


1.718 


71 


130.3806 


7.449 X 


10-189 


1.713 


72 


132.0878 


3.173 X 


10-185 


1.707 


73 


133.7906 


1.306 X 


10-181 


1.702 


74 


135.4867 


5.073 X 


10-178 


1.696 


75 


137.1783 


1.906 X 


10-174 


1.691 


76 


138.8633 


6.765 X 


10-1^1 


1.684 


77 


140.5439 


2.322 X 


10-167 


1.680 


78 


142.2178 


7.550 X 


10-164 


1.673 


79 


143.8850 


2.345 X 


10-160 


1.667 


80 


145.5478 


6.996 X 


10-157 


1.662 



No. 


Energy 


Penetrability, 


AE 




Erad 


Tb 


ar 




81 


147.2040 


1.990 X 


10-15^5 


1.656 


82 


148.8534 


5.400 X 


10-150 


1.649 


83 


150.4985 


1.409 X 


10-146 


1.645 


84 


152.1368 


3.518 X 


10-143 


1.638 


85 


153.7686 


8.353 X 


10-140 


1.631 


86 


155.3936 


1.900 X 


10-136 


1.625 


87 


157.0120 


4.109 X 


10-133 


1.618 


88 


158.6260 


8.619 X 


10-130 


1.613 


89 


160.2310 


1.700 X 


10-126 


1.605 


90 


161.8295 


3.224 X 


10-123 


1.598 


91 


163.4234 


5.873 X 


10-120 


1.593 


92 


165.0085 


1.018 X 


10-116 


1.585 


93 


166.5870 


1.679 X 


10-113 


1.578 


94 


168.1587 


2.666 X 


10-110 


1.571 


95 


169.7216 


3.987 X 


10-107 


1.562 


96 


171.2801 


5.801 X 


10-104 


1.558 


97 


172.8274 


7.876 X 


10-101 


1.547 


98 


174.3704 


1.040 X 


io-»7 


1.542 


99 


175.9044 


1.297 X 


10-94 


1.534 


100 


177.4318 


1.557 X 


io-«i 


1.527 


101 


178.9480 


1.750 X 


10-88 


1.516 


102 


180.4599 


1.911 X 


10-85 


1.511 


103 


181.9606 


1.957 X 


10-82 


1.500 


104 


183.4547 


1.925 X 


10-^9 


1.494 


105 


184.9377 


1.782 X 


10-^6 


1.482 


106 


186.4140 


1.581 X 


lo-^-'^ 


1.476 


107 


187.8814 


1.337 X 


10-^" 


1.467 


108 


189.3378 


1.060 X 


10-67 


1.456 


109 


190.7874 


8.123 X 


10-65 


1.449 


110 


192.2238 


5.763 X 


10-62 


1.436 


111 


193.6535 


3.957 X 


10-59 


1.429 


112 


195.0699 


2.518 X 


10-56 


1.416 


113 


196.4773 


1.534 X 


10-53 


1.407 


114 


197.8737 


8.777 X 


10-51 


1.396 


115 


199.2590 


4.742 X 


10-48 


1.385 


116 


200.6332 


2.420 X 


10-45 


1.374 


117 


201.9963 


1.160 X 


10-42 


1.363 



Table 3: The resonant levels inside the energy region -E^ad = - 200 {AE = -Erad.noxt — Srad.prov, we have choose: 
csstart ~ 10, Omax = 70, number of intervals inside the scale axis a is 1000, number of intervals of energy is 100000). 
We also add the coefficients of penetrability Tbar, oscillations Kosc and penetration Tmir for such levels (we used 
A = 0.001, B — 0.001). For all presented energy levels, -Erad, we have achieved accuracy in calculations of Tbar and i?bar: 

|Tbar+i?bar-l| < I-IQ-^^ 
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